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Multi-high-frequency electron paramagnetic resonance (EPR) spectrum 
for a super molecular dimer [Mn^ of single-molecule magnets recently re- 
ported [S. Hill, R. S. Edwards, N. Aliaga-Alcalde and G. Christou(HEAC), 
Science 302, 1015 (2003)] is studied in terms of the perturbation method in 
which the high-order corrections to the level splittings of degenerate states are 
included. It is shown that the corresponding eigenvectors are composed of en- 
tangled states of two molecules. The EPR-peak positions are calculated in 
terms of the eigenstates at various frequencies. From the best fit of theoretical 
level splittings with the measured values we obtain the anisotropy constant 
and exchange coupling which are in agreement with the corresponding val- 
ues of experimental observation. Our study confirms the prediction of HEAC 
that the two Mn^ units within the dimer are coupled quantum mechanically 
by the antiferromagnetic exchange interaction and the super molecular dimer 
behaviors in analogy with artificially fabricated quantum dots. 
PACS numbers: 75.45.+j, 75.50.Xx, 75.50.Tt 



I. INTRODUCTION 



Single-molecule magnet s(SMM), which may be the smallest nanomagnets exhibiting mag- 
netization hysteresis- loop ( a classical property of macroscopic magnets), straddle the inter- 
face between the characteristics of classical and quantum worlds. The quantum tunneling of 
magnetization and quantum phase interference known as macroscopic quantum effects have 
become an attractive research field in recent years. The quantum effects in these magnetic 
particles and clusters of nanometer size, such as Mnu(S = 10),Fe 8 (5' = 10) and Mn^S = |) 
molecules, have been well studied 1-17 both experimentally and theoretically. Several propos- 
als of possible quantum computing schemes are also suggested using molecular magnets 18-20 . 
The high-frequency single crystal electron paramagnetic resonance (EPR) is a powerful tool 
for the experimental study of various properties of single-molecule magnets 21-24 . Werns- 
dorfer et al. 25 recently pointed out that the supermolecular dimer [Mr2, 4 ] 2 consisting of two 
molecule magnets Mn 4 with antiferromagnetic exchange-coupling exhibits a quite different 
quantum behavior from two individual Mn 4 molecules without coupling 25-28 . It is there- 
fore of great importance to understand the effect of the exchange interaction which leads to 
energy level splitting 27 ' 31 and has been studied with the high-order perturbation method 29 . 
The step-like magnetization hysteresis-loop in a supermolecular dimer have been demon- 
strated in terms of the numerical solution of the time- dependent Schrodinger equation 30 . 
Recently the multi-high-frequency EPR was used to probe the magnetic excitations of the 
supermolecular dimer [Mn 4 ] 2 by Hill et al 31 . The measured spectra well interpret both the 
quantum transitions involving coherent superposition states of the two molecules and the 
phase decoherence rate, which provide a compelling evidence that the molecules are coupled 
quantum mechanically by the antiferromagnetic exchange interaction. In this paper, we 
restudy the EPR transitions in the supermolecular dimer [Mn4] 2 employing the high-order 
perturbation method 29,32 and obtain the high-order corrections of level splitting. In terms 
of the corresponding eigenvectors accurate values of the anisotropy constant D and the an- 
tiferromagnetic exchange-coupling J are determined by the best fit between the theoretical 
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and measured level-splitting. The EPR- peak positions are also calculated as a function of 
frequencies. 

We in the following section first construct the explicit eigenvectors which are seen to 
be composed of entangled states. It is then shown that the ground states of the dimer are 
nearly (with probability 0.99883) the maximum entangled states, or Bell states which are 
useful in the quantum computing. The main goal of this work is to display all eigenvalues 
and the entanglements of the corresponding eigenstates in the dimer system. Our theoretical 
studies shine more light on the prediction of Hill et al 31 that the supermolecular dimer with 
proper chemical design may be used as practical quantum devices for quantum computing. 



II. LEVEL SPLITTING AND EIGENVECTORS 

Neglecting off-diagonal crystal field terms and inter-molecular interactions, the effective 
spin Hamiltonian with a magnetic field applied parallel to the easy axis of an isolate SMM 
has the form 31 Hi = —D Sf z — B^O® + g z n B B z S iz , .where S iz is the z-axis spin projection 
operator, and the index i (= 1,2) is used to label the two Mn 4 molecules in the dimer 
; D > denotes the uniaxial anisotropy constant; B^O® characterizes the fourth order 
uniaxial anisotropy where 0° = - [30S(S + 1) - 25]S 2 Z - 6S(S + 1) + 3S 2 (S + l) 2 ; g is 
the electronic g-factor; and hb is the Bohr magneton. The Hamiltonian Hi may be written 
in the following form apart from a trivial constant 

Hi = -DSl - BSf z + g z fi B B z S iz , (1) 

where D = D — 30S(S + 1) - 25, B = 355°. 

For the supermolecular dimer with antiferromagnetic coupling the Hamiltonian is seen 
to be 



H 



H ± + H 2 + J z S lz S 2z + -J xy (S 1+ S 2 - + 5'i_S i 2 + ) , (2) 

where S± — S x ± i S y with J z and J xy denoting the strengths of exchange coupling. In the 
following we consider only the isotropic case that J z = J xy = J > which has been verified 



experimentally 27 . We furthermore assume that the exchange coupling J is much less than 
the anisotropy constant D. Then, H is rewritten as the perturbation form 

H = H + V, (3) 
H = -D (S 2 lz + Si) - B (Sf z + Si) + JS lz S 2z + g gflB H g (§ lz + S 2z ) , (4) 

V = l -J (S 1+ S 2 _ + Si_S 2+ ) , (5) 

where Hq is the zeroth-order Hamiltonian. The zeroth-order eigenvectors of the dimer are 
direct products of the single-molecule eigenvectors such that \S 1 ,S 2 ,M 1 , M 2 ) or abbreviated 
as|Mi,M 2 ), where Mi and M 2 represent the quantum numbers of spin operators S\ z and 
S 2z respectively. The dimer has (2 Si + 1)(2S 2 + 1) energy levels which are labeled by the 
spin quantum numbers Mi and M 2 and given by 

E Q = -D (Ml + Ml) - B (M? + M 2 4 ) + JMiM 2 + g zl i B H z {Mi + M 2 ) . (6) 

It is seen obviously that the eigenstates \Mi,M 2 ) and \M 2 ,Mi) are degenerate since 
Eq(Mi, M 2 ) = Eq(M 2 , Mi). The two Mn\ molecules in the [Mn^ dimer are coupled by a 
weak exchange interaction via both the six C-H- • -CI hydrogen bonds. Thus, we can treat 
V perturbationally.. The level splitting of the pair of degenerate eigenstates \Mi,M 2 ) and 
\M 2 , Mi) with M 2 > Mi appears only in the chain of matrix elements connecting the states 
| Mi + k, M 2 - k) and \M ± + k + 1, M 2 - k - 1) ,where k = 1, 2, • • • , M 2 - Mi - 1. Then, the 
level splitting becomes 29 ' 32 

AE MlM2iM2Ml = 2Vm 1 m 2 ,(Mi+i)(m 2 -i)t; f; V(m 1 +i)(m 2 -i),(m 1 +2)(m 2 -2) 

^(M 1 +l)(M 2 -l) ~ &M!M 2 
p V r (Af 2 -l)(Afi+l),M- 2 Mi ) (7) 



E(M!+2)(M 2 -2) — Em x M 2 



where 



Vm 1 m 2 ,(m 1 +i)(m 2 -i) = (Mi,M 2 | - J (§i + S 2 - + Si-S 2+ ) \Mi + 1,M 2 — 1) = — h( Ml +i)(M 2 -i), 

(8) 



and 



h(M 1 +l)(M 2 -l) — [(Si + Mi + 1) (Si - Mi) (S 2 -M 2 + l) (S 2 + M 2 )Y 



(9) 



Since Si — S 2 — S for the dimer case, we obtain the level splitting as 

w jMi-m, (S + M 2 )\ {S - Mi)\ 



AE = \ II [ 2D ' - ABk ( M 2 — Mi — k) 



[(Mi - M 2 - l)!] 2 (5 - M 2 )! (5 + Mi)! 



, (10) 



where 



D' = 2D + J + AB (Mi + M 2 f - ABMiM 2 . 



(11) 



Disregarding the fourth- power of anisotropy term in Eq..(4) i.e. B = 0, the formula of level 
splitting can be simplified as 



A£ = < 4fl + 2J »(«JT2j) 



M2-M1 



(S + M 2 )\(S-Mi)\ 



1 



(M 2 - Mi - 1)! 



(12) 



AD + 2JJ (S - M 2 )\ (S + Mi)\ 

The state of the zeroth-order perturbation is obviously the superposition of the states 
|Mi,M 2 ) and |M 2 ,Mi) such that 



^ (0) ) = [Mi, M 2 ) + a { S 2Ml \M 2 , M ± ) . 



(0) 



(13) 



The first-order state is then obtained as 



(i) 



E a kik2 \ k i,h), 

ki,k 2 



where 



(0) T> 1 (0) fr 

(1) _ a M 1 M 2 v k 1 k 2 M\M 2 + a M 2 M 1 Vk\k 2 MiM\ 



a kik 2 ~ 



E 



M 2 M\ — E^k 2 



and OmiM 2 = ^Imi = 0- The n-th order state is 



where 



ki,k 2 



(n) _ \ - ^fcifc 2 ,fc^ 
a fcifc 2 — P P Vfe' ' 

,/ ,/ &m 2 M x - ^k x k 2 12 

1 ' 2 



(14) 



(15) 



(16) 



(17) 
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and (ki, fc 2 ) ^ (fci, k' 2 ) ^ (M l ,M 2 ) ^ (M 2 , Mi) , afi^ = ag Mi = 0(m = 1, 2, • • • , n). For n 
= M 2 — k\ — k 2 — Mi > 1, we have 



» 



^fcifc 2 ,(fci+l)(fc2-l) ^(fci+l)(fc2-l),(fci+2)(fc2-2) V(fci+2)(fc2-2),(fci+3)(fc 2 -3) V(M 2 -l)(Mi+l), 



a kik 2 i- i i- i- 

&M 2 M 1 — ftkiki Hi Mi Mi ~ 



M 2 M\ 



J k\k 2 



fel+l,fc 2 -l 



E 



M 2 Mx 



E, 



fci+2,fc 2 -2 



E 



M 2 M\ 



Em 2 -i,m 1 +i 
(18) 



and 



J 



^(fci+/-i)(fc2-i+i),(ifci+0(fc2-0 — -^{ki + l — l,k2 — l + l\ Si-S 2 + \h + l,k 2 — 1} 

= J ~ [(S - k 2 + (S + k 2 + 1 - (S + h + (S - h + 1 - (19) 

where I = 1, 2, • • • , n. We finally obtain the coefficients of the n-th order state 



tX = r { II t M i2 + Z (A;i - fc 2 + 0] f2D" + ABl (k, -k 2 + l) 



'n—l 



-1 



. 1=0 



X 



(5 + M 2 )\ (S - M l )\ {S - h)\ (S + k 2 )\ 



(20) 



(S - M 2 )\ (S + Mi)! (5 + fcO! (5 - A; 2 )!_ 

where M i2 = M X M 2 - k x k 2 , D" = 2D + J + AB (Mi + M 2 ) 2 - 2S (M X M 2 + j^fe) . For the 
case n = k[ — Mi = M 2 — k' 2 > 1 ,we have 



n -1 



^ = J n ' J J] [M{ 2 + / (fc 2 - fci + /)] [2D'" + 45/ (k 2 -k[ + l 



x 



(5 + M 2 )!(5-M!)! (5-4)! (S + k[)\ 



(21) 



(5-M 2 )!(S , + Mi)!(5 + *4)!(S'-A;i)! 

where M{ 2 = M X M 2 - k[k' 2 , £>"' = 2D + J + 4B (Mi + M 2 ) 2 - 25 (m x M 2 + fc^) . 

The zeroth-order perturbation wave functions for the non-degenerate states |Mi,M 2 ) 
with M M> M are seen to be ^ (0) ^> = 2 |M, M) . The matrix elements of op- 
erator V in Eq. 18 appear only in the chain connecting the states \k± + l,k 2 — I) and 
\k± + I + 1, k 2 — I — 1) . We can obtain the results by non-degenerate perturbation theory 



as follows. When a^ k2 



+ 0,a 



(i) 



,( 2 ) 



feifc 2 ~~ a k\k 2 



(n-i) n 

= a klk , 2 = 0, as a consequence eqs. 



(16), (17), (20) and (21) become the same as for the case of non-degenerate states. 
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Making use of eqs.(13), (16), (20) and (21), we can obtain all eigenvectors for the quantum 
number range from M = Mi+M 2 = —9 to , while the eigenvectors for the quantum number 
range from M — 1 to 9 can be obtained simply by the replacements: M\ — > —Mi,M 2 — > 
— M 2 in the former case. The eigenstates are listed in appendix A (where \M 1 ,M 2 ) S and 
\M 1 ,M 2 ) A denote the symmetric and antisymmetric states respectively ) with parameter 
values D = 0.72K, B = 1.8 x 10~ 3 K, J = 0.10-fT which give rise to the best fit between 
theoretical and measured level splitting. The parameters D and B determined here with 
the best fit are in agreement with the values obtained by Hill et al 31 . When Mi ^ M 2 , 
the normalized symmetric \Mi,M 2 ) s and the antisymmetric \Mi,M 2 ) A states are actually 
maximum entangled states with the entanglement degree of von Neumann entropy that 
E(p) = 1, where 



and p is the reduced density matrix 33 . When Mi = M 2 , the states are disentangled. Thus, 
the eigenvectors are, in general, composed of maximum entangled states with definite prob- 
abilities. For example, \ifii) s is the disentangled state &nd\ip 2 ) s , \ip2) a > \^)a are maximum 
entangled states with E(p) = 1. The state \ip^) A is composed of maximum entangled states 



states 1^26)5^4 is 0.99883 . It is also seen that the ground states ^26)5 and \^2%) A are vei T 
close to the maximum entangled states in a small range of magnetic field values. 



In the EPR experiment, it has been seen that the exchange-coupling induced single- 
spin transitions i.e. the transitions from (M l5 M 2 ) to (M x + 1,M 2 ) or from (Mi, M 2 ) to 
(Ml,M 2 + 1) can occur, however, the single-spin transitions should also depend on the 
states of the other spin within the dimer because of the exchange coupling. This depen- 
dence has displayed its importance in the resonant macroscopic quantum tunneling too. 



E(p) = Trplog 2 p, 




with probabilities 0.9445 and 0.0555 respectively. More- 



Ill. EPR TRANSITION 
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Most importantly the multi-high-frequency electron paramagnetic resonance spectra of the 
dimer obtained by Wernsdorfer et al 27 and Hill et al. 31 show a compelling evidence that the 
two molecules Mn 4 are coupled quantum mechanically by the antiferromagnetic exchange- 
coupling which leads to the energy level splitting.. In terms of Eq.(6) and Eq. (10), we 
can calculate the energy level splitting caused by the perturbation part of the Hamiltonian 
V. Thus, when the frequency v is fixed to the value lAhGHz and the magnetic field B z 
varies quasi-statically to reach a resonance point, the EPR transition matrix elements with 
the selection rule AM = ±1 can be accurately calculated using the level splitting formula 
Eq. (10) with the help of the corresponding eigenvectors shown in appendix A. The mag- 
netic dipole perturbation allows transitions only between states of the same symmetry. The 
positions of principal and other possible EPR peaks are given in Table 1 (resonance tran- 
sitions between higher energy levels are not listed ), and displayed schematically in Fig. 1. 
From Table 1, we can find that principal resonance peaks labeled by (x), (a) to (i) are in 
agreement with the result of Hill et al. 31 (see Fig. 2 and Fig. 3 in Ref. [31]). From the 
Fig. 3 of Ref. [31] we can see that the resonance peaks labeled by (/), (m), (n), (p) and (q) 
with magnetic field values B z = 2.75 ± 0.05, 4.85 ± 0.05, 2.20 ± 0.05, 3.55 ± 0.05 , 3.80 ± 0.05 
respectively, are also in agreement with our result in Table 1. We expect that the possible 
EPR peaks calculated at the higher field region for B z > 6 tesla can be verified by future 
EPR experiments. Moreover the relation between the field value B z and frequency v can 
be obtained from the resonance condition 

hv= E(M + 1)- E(M) = AE + g z fi B B z , (22) 

where the energy E (M) can be accurately calculated from the level splitting formula Eq. 
(10) and the eigenvectors in appendix A. The level splitting AE is shown in Table 2 for 
various EPR transitions. The calculated positions of EPR peaks depending on frequencies 
are shown in Fig. 2. 
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IV. CONCLUSION 



The EPR spectra in the supermolecular dimer [Mn 4 ] 2 with the antiferromagnetic 
exchange-coupling between two molecule magnets are studied in terms of the high-order 
perturbation method which allows us to obtain the level splitting of the degenerate pairing 
states and the corresponding eigenvectors. It is shown that the eigenvectors are composed 
of maximum entangled states with definite probabilities. This observation shines more light 
on the potential application of molecular magnets in quantum information. The theoretical 
results are in good agreement with the measured values 31 . The EPR peaks are obtained 
in terms of the high-order level splitting formula and eigenvectors from which the value of 
exchange coupling J = 0.10K is determined and is seen in accord with the result of the reso- 
nance quantum tunneling 25-29 . Both experimental and theoretical investigations of EPR for 
the supermolecular dimers are of great importance for a better understanding the dynamics 
of such systems which are potential devices for quantum computing. 
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Appendix A: Eigenvectors 
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Table caption 

Table 1. The resonant field transitions. Theoretical values are obtained by the high-order 
level splitting formula and eigenvectors with parameter values D = 0.72K, B = 1.8 x 10~ 3 K 
and J = 0.10K. 

Table 2. Level splitting AE in Eq. (22). 

Figure caption 

Fig. 1. Schematic display of lowest energy states (M= —9 to —4) and the part of lower 
energy states (M= —3 to 1). The zeroth-order energy levels and eigenvectors are shown on 
the left; energy shifts due to the exchange-interaction are shown in the center; and the results 
of the high-order perturbation calculation are displayed on the right. The red lines [x, (a) to 
(i)] denote the strongest EPR transitions; the blue lines denote the weaker EPR transitions; 
the magenta lines denote other possible EPR transitions. 

Fig. 2. EPR- peak positions ( solid lines) obtained from Eq. (23) with different frequen- 
cies. 
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